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ABSTRACT 

p/| We report the discovery of an extended halo (~ 40' in diameter) around the 

C/3 ■ planetary nebula NGC 7293 (the Helix Nebula) observed in 12 jum band from 

the Wide-field Infrared Survey Explorer all-sky survey. The mid-infrared halo 
has an axisymmetric structure with a sharp boundary to the northeast and a 
more diffuse boundary to the southwest, suggesting an interaction between the 
stellar wind and the interstellar medium (ISM). The symmetry axis of the halo 
is well aligned with that of a northeast arc, suggesting that the two structures 
^ . are physically associated. We have attempted to fit the observed geometry with 

a model of a moving steady-state stellar wind interacting with the ISM. Possible 
combinations of the ISM density and the stellar velocity are derived from these 

SO 



Of 



O 



fittings. The discrepancies between the model and the observations suggest that 
the stellar mass loss has a more complicated history, including possible time and 



angle dependences. 

Subject headings: infrared: ISM — ISM: structure — planetary nebulae: indi- 
vidual: Helix Nebula (NGC 7293) — stars: AGB and post-AGB 

1. INTRODUCTION 

As low- and intermediate-mass stars evolve from asymptotic giant branch (AGB) to 
planetary nebula (PN) stages, they continually eject matter into the interstellar medium 
(ISM). The resultant circumstellar envelopes are one of the main sources of chemical en- 
richment of the Galaxy. Our current understand ing of the origin of PNs is based on the 



interacting stellar wind (ISW) model ( lKwoklll982l ) which has predicted the presence of AGB 



halos created by AGB winds. The remnants of this AGB wind are swept up by a later- 
developed fast wind to form a PN. One of the predictions of the model is that the remnants 
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of the AGB wind should be detectable as halos in PNs. Since then, AGB halos have been 
detected by a variety of techniques. These include observations (1) in the visible through 
emission or scattering, (2) in millimeter and submillimeter wavelengths through molecular- 
line emission, and (3) in the infrared through dust emission. Since PN halos can contain more 
mass than the main bright nebular shells, studies of PN halos are important to understand 
the AGB mass-loss history and the chemical evolution of the Galaxy. 

In the visible, AGB halos are usually very faint (with surface brightness about 10~ 3 - 
10~ 4 lower than the main nebulae), making th em difficult to observe. Many efforts have been 



made to search for faint AGB halos (see e.g. ICorradi et al.ll2003l ; iRamos-Larios &: Phillips 



20091 ). These observations show that the detected AGB halos usually exhibit a bow-shaped 
bright rim, suggesting the interaction between AGB winds and the ISM. The wi nd-ISM 



inter a ction process has been the subject of various theoretical stu d ies (see e.g. ISmith 



1976 



Borkows 



Wareing et al. 



ri et a" 



2007a 



199p[ ISoker et all Il99ll : bgani & Soke] Il998t Ivillaver et all I2OO3I : 



b|). Early studies suggest that the wind-ISM interaction is observ- 



able only for evolved PNs whose density has dropped below a critical limit. More recent 
hydrodynamic simulations argue that the interaction can be observed for young PNs if the 
time-dependent variations in AGB and post-AGB winds are t aken into account. The w ind- 
ISM interaction was indeed observed in not only a ncient PNs (ITweedv K witter! 19961) but 



also in the circumstellar envelopes of AGB stars ( lUeta et al. 
2012h . 



2006 



Ueta 



2008 



Cox et al. 



The high-quality images acquired by recent wide field-of-view (FOV) surveys have been 
pro ven useful to search for l arge-s cale extended structures around evolved stars. For exam- 
ple, ISahai fc Chronopoulosl ( 120101 ) discovered a large outer shell around the highly evolved 
AGB star IRC+10 216 using the images taken with Galaxy Evolution Explorer (GALEX; 
Martin et al.l 120051 ) . In this paper, we report the results of searching for outer structures 
around NGC 7293 utilizing the recently released m id-infrared images of the Wide-field In- 
frared Survey Explorer ( WISE; [Wright et al.ll2010f ) all-sky survey. 



NGC 7293 is one of the most extensively s tudied PNs (see lO'Dell et al.ll2004j ; iHora et al. 



20061 : ISpeck et al.l 120021 : Meaburn et al.lll992l I2005L I2008L and reference therein). Given its 
relative proximity to the Sun and its relatively accurate distance measurement, NGC 7293 
is an excellent target to study in detail the spatial structu re of PNs. In th i s pap er, we have 



adopted the distance estimate of 216_ 12 pc as determined by lBenedict et al.l ( 120091 ) . Although 



its overall appearance seems simple, a number of intriguing structures, from small cometary 
knots to large-scale arcs, to bipolar outflows, have been revealed in previous observations. 
The inner helical structure is composed of thousands of cometary knots of lowly-ionized and 
molecular gas. The main nebula consists of two rings of highly-ionized gas and a faint outer 
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fil ament. The thr e e-dim ensio nal (3-D) structure of the main nebula has been investigated 
by iMeaburn et al.l (120081 ) and lO'Dell et al.l (12004 ) who suggested that it was formed through 
multiple outflows along different directions. 

Although NGC7293 is of great astrophysical interest, it is not an easy object to ob- 
serve because of its low surface brightness and its large angular extent in the sky. Without 
highly sensitive and large FOV observations, it is impossible to explore the faint halo around 
this PN. The recent WISE observations, therefore, provide an invaluable resource for that 
purpose. 



2. The Data 

The data used in this paper were taken from the WISE All-Sky Data Releas^E The data 
have been calibrated with the VF/ST^-pipeline. WISE mapped the entire sky in four bands 
centered at wavelengths of 3.4, 4.6, 12, and 22 /im (W1-W4) with an angular resolution of 
6.1", 6.4", 6.5", and 12.0", respectively. During the survey, the four bands were simultane- 
ously deployed to image a 47' x 47' FOV on the sky. The exposure time of individual frame 
is 7.7 seconds in Wl and W2, and 8.8 seconds in W3 and W4. The number of exposures are 
different at different ecliptic latitudes. The observations of the fields including NGC 7293 
were made on 20-23 May 2010, with total accumulated time of 17.2 minutes for Wl and 
W2, and 19.7 minutes for W3 and W4. The composite-color WISE image of NGC 7293 is 
shown in Fig. [TJ 

Our analysis is supplemented by the far- and near-ultraviolet (FUV and NUV) images 
retrieved from the GALEX archive. The FUV and NUV images have effective wavelengths 
of 1516 A and 2267 A, and angular resolutions of 4.5" and 6", respectively. A detailed de- 



script ion of the calibration and data products of GALEX can be found in IMorrissey et al. 



( 120071 ). The observations of NGC 7293 were made on 22-25 August 2004 as part of the 
Calibration Imaging Survey (CAI). The accumulated exposure time is 56 minutes for each 
band. In order to improve the signal-to-noise ratio, we co-added the FUV and NUV images. 
The resultant GALEX image is shown in the right panel of Fig. [TJ 



http://wise2.ipac.caltech.edu/docs/release/allsky/ 
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3. RESULTS 

An axisymmetric extended halo at 12 /mi surrounding NGC 7293 can clearly be seen in 
Fig. [TJ The halo exhibits a bow-like shape with a well-defined boundary to the northeast 
(NE), in contrast to a diffuse and open structures to the southwest (SW). Such asymmetric 
morphology gives a clear indication that the object is moving toward the NE. The direction of 
motion can be approximated by the symmetry axis of the halo and the symmetry axis is found 
to located at a position angle (PA) of 60°. The NE boundary of the halo can be approximately 
described by a circle with a radius of 1213" (3.9 x 10 13 km= 2.6 x 10 5 AU= 1.3 pc), with a 
center located at 359" (1.2 x 10 13 km= 7.8 x 10 4 AU= 0.4 pc) from the central star (see Fig. 

We find that the 12 /im emission in the halo is relatively smooth and is about 10 times 
fainter than that of the main nebula. This halo is not detectable at the WISE 22 /mi band 
and has not been seen in previously reported optical images^. It is only marginally detectable 
in the WISE 3.4 and 4.6 /im and in the GALEX images. Assuming that the infrared emission 
peaks at 12 /im, we infer that the emission of the halo is dominated by thermal emission 
from warm dust with a temperature of ~ 300 K. One can in principle obtain the spatial 
distribution of temperature in the halo by comparing the W1-W3 images. However, the 
halos in Wl and W2 are too faint to yield meaningful conclusions. Nevertheless, we cannot 
rule out the possibility that shock-excited H 2 li nes might significant l y cont ribute to the mid- 



infrared emission of the halo, as suggested by iRamos-Larios et al.l ( 120101 ) in their study of 
the halo of NGC 40. 

The shape of the halo obviously suggests that it is interacting with the IS M. The conclu- 



sion is strengthened by the study of the proper motion of the Helix Nebula by iKerber et al. 



(120081 ). The direction of the proper motion is found to be at PA= 77°, which roughly aligns 
with the symmetry axis of the halo and the bow apex. The apex of the bow-shaped halo 
(position A in Fig. [2]) is about 854" (2.8 x 10 13 km=1.8 x 10 5 AU=0.9 pc) from the central 
star. Assuming the velocity of AGB wind, v w = lOkrns -1 , we obtain a lower limit to the 
dynamic age of the halo of 88 000 yr, suggesting that the halo was inde ed formed long ago , 



and can therefore serve as a fossil record of the AGB mass-loss history. lO'Dell et al.l (120041 ) 
estimated that the outer material of NGC 7293 has a recombination timescale of ~ 40 000 yr, 
much shorter than the dynamic age of the halo. Therefore, it is very likely that the halo is 



neutral and optically invisible. Based on optical observations, iMeaburn et al.l ( 120081 ) deter- 



mined the dynamic age of the main nebula to be llOOOyr. If the halo has been continually 



2 After this paper was published, we were informed by Q. Parker that their deep Ha image has revealed 
some faint extended structures that are associated with the halo (Parker et al. 2001). 



-5 - 



developed with a mass-loss rate of M — 10~ 5 M yr _1 and was terminated when a fast wind 
was launched to form the main nebula, the mass of the halo can be estimated to be >O.77M . 
Due to the wind-ISM interaction, the past of the halo material has probably been stripped 
to the SW behind the central star. 

The WISE images also reveal some smaller-scale structures, such as the NE arc and 
small cometary knots. The main nebula exhibits a collimated bipolar outflow oriented NW- 
SE (the NW and SE outer features in Fig. 1, gned by O'Dell et al. 2004). We 

note that the SE ends of the bipolar outflow are running into the boundary of the halo. 
These results are consistent with the n a rrowe r FOV infrared observations acquired with the 
Spitzer Space Telescope by iHora et al.l (120061 ). The NE arc is the brightest feature outside 
the main nebula. Fig. [2] shows that the 12 fim image of the NE arc can be described by a 
circle arc about 696" from the central star while its FUV emission is closer to the central 



star (about 680" away). IHora et al.l (120061 ) found that the H2 mid-infrared emission of the 
NE arc is located in a shell outside the Ha emission. Because the GALEX FUV emission 
is presumably from shock-excited H2, the stratified optical/FUV and mid-infrared emission 
structures likely suggest the radial variations of excitation conditions rather than the H + /H 2 
ratios. We also note that the WISE 22 /zm image shows strong brighte ning around the 
central star. This confirms the previous Spitzer observations at 24 /zm by ISu et al.l (120071 ) 
who ascribed the 24 /im brightening to the presence of a inner debris disk. 

In the outer region, the GALEX and WISE images appear strikingly different (see Fig. 
H]). The GALEX image shows many bow-shaped features that are not dete cted by WISE. 



A sch ematic presentation of the observed structures is presented in Fig. 6 of iMeaburn et al. 



(120081 ). These FUV features are presumably due to H 2 emission excited by fast shock waves. 
The most interesting feature is the NE bow-shaped filament, which contains three bow- 
shaped sub-structures (61-63 in the inset of Fig.[T]) and whose apex is exactly aligned with the 
direction of proper motion (see the proper motion vector in Fig.[T]). These sub-structures were 
pr obably formed due to instabilities in the shock interface, as predicted by the simulations 
of iMohamed et al.l (120121 ). We note that the apex of the fast outflow shown in the GALEX 
image has a "L" shape, which is characteristic of an outflow near the plane of the sky. If the 
outflow is inclined at a large angle, the apex would take on a more rounded shape. 

Unlike the bow-shaped filament, the NE boundary of the halo cannot be seen in the 
GALEX image. This is probably the result of the shock responsible for the halo boundary 
being too slow to provide enough e nergy to excite Ho emission. The GALEX image also 
shows the NE object and shock (see lO'Dell et al.ll2004 for a discussion), as well as bipolar 
streams in NE-SW orientation (marked as "outer features" in Fig. [TJ). Although there are 
several asymmetrical outflows (outside and inside the halo), the halo boundary seems to 
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be continuous and not broken (as seen in the left panel of Fig. 1). This simply suggests 
that from the perspective perpendicular to the line of sight, these streams do not cross the 
projected halo boundary. 

The mid-infrared halo has a appearance similar to the shell around IRC+10216 detected 
by ISahai fc Chronopoulosl (120101 ) . except that the vortices behind the IRC+10216 shell are 
not clearly detected in NGC 7293. If the vortices had been formed behind the halo of 
NGC 7293, they supposedly would have been diluted to the ISM. 

Overall, the newly discovered mid-infrared halo of the well studied PN NGC 7293 
strengthens the conclusions that AGB halos are commonly pres ent and infrared imaging 



is pa r ticularly useful to detect the faint extended structures (e.g. iRamos-Larios Phillips 
20091 ; IRamos-Larios et al.l 120121 ) . The direct detection of the halo of NGC 7293 provides us 
with the knowledge of the ambient environ ment within which the later-developed dynamical 
processes [e.g., the multipolar fast outflows ( jMeaburn et al.ll20081 )] operates. By deriving the 
geometry and density of the halo, we are in a much better position to constrain the physical 
processes at work. 



4. DISCUSSION 



The shape of the 12 /xm halo suggests the presence of a bow shock. Since the shape 
of bow shock is dependent on the velocity and density contrasts between the matter of 
the stellar wind and the ISM, the observed structure can serve as a useful probe of the 
properties of the AGB wind and the ISM. If the post-shock material has a high cooling rate, 
the physical thickness of the shell at the shock interface would be small. Under the thin-shell 
approximation and assuming t hat an equilibr ium between the ram pressure of stellar wind 
and that of the ISM is reached, IWilkinl ( 119961 ) deduced an analytic expression for the shape 
of the bow shock, 



R(6) = Rq sin" 1 0^/3(1 -0 cot 0), (1) 

where is the polar angle from the symmetric axis as seen from the central star at the 

coordinate origin, and Rq is the distance between the apex and the central star (i.e. stand- 
off distance) given by 

Ro = \jMv w /{Airp lSM vl). (2) 

In eq. 2, pism is the mass density of the ISM and the velocity of star respect to the ISM. 
R{9) increases with increasing because the ram pressure of the ISM is the strongest at 
the apex and decreases on either side. Assuming that the properties of the stellar wind are 
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typical (see Tabled]), one can obtain the relation between v* and the number density (tiism) of 
the ISM from eq. (2) by measuring Rq. Wilkin's analytic solutions have been widely applied 



supergiants, AGB stars, and pulsars (e.g. 
201 lh . 



Cox et al. 


2012; 


Mackev et al. 


2012; 


Jorissen et al. 



We have attempted to determine the R value by fitting the projected shape of a 3- 
D paraboloid (constructed from eq. 1) on the sky plane with the observations. For given 
assumed values of Rq and i (inclination angle of the symmetric axis respect to the sky plane; 
no matter whether the apex points toward or away from us), the projected distances between 
every point on the paraboloid and the central star can be determined. The boundary of the 
projected 3-D paraboloid is found by calculating the furthermost projected distance at each 
position angle in the sky plane. Valu es of Rp and % ar e then adjusted to obtain the best fit. 
Similar procedures have been used by lCox et al.l (120121 ) to study the bow shocks around AGB 
stars and red supergiants. We should note that since the bow shock has a sharp boundary, 
the inclination angle cannot be too large (i.e., the star is moving not far from the plane of 
the sky). 

We constructed a series of models with different % and Rq values, but found none of them 
can completely account for the shape of the halo boundary. Fig. [2] shows four of the models 
and the model parameters are listed in Table[TJ We can see that as i increases, the projected 
location of the apex of the paraboloid changes from A toward the central star (positions a 
and b). If i were 90°, the projected position of the apex would be at the central star. The 
models suggest that given a certain distance from the projected apex (position A) to the 
central star, the values of Rq monotoni cally decreases with in crea sing i (see model s 1 and 
2). This is consistent with the results of iMac Low et al.l ( 119911 ) and I Cox et al.l (120121 ) . Since 
the shapes of the model curves are insensitive to i, we are unable to match the observed 
positions of B and C. If we force the model curves to pass through B and C (see models 3 
and 4), then the model boundary curves deviate from the position of A. Nevertheless, our 
models suggest that the t>*/ y'nisM value is probably in the range of 15-50 km s _1 cm 1 - 5 . For 
example, in model 1, the number density of H in the ISM would be 0.3 cm" 3 if the star is 
moving at a speed of 10 km s _1 . 

Below we will discuss the possible causes for the inability of Wilkin's expression in fitting 
the observations. If the post-shock material is adiabatic, a termination shock and a bow shock 
are expected to form inside and outside the halo, respectively. In this case, the boundary of 
the halo is the discontinuity between stellar and interstell ar material (i.e. astropause). These 
structures have been detected in the shell of IRC+10216 (jSahai fc Chronopoulosll2010l ). One 
may have an impression from Fig. [1] that the NE bow-shaped filament detected by GALEX 
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and the NE arc are located outside and inside the halo respectively, and perhaps they 
represent the bow shock and the termination shock. This scenario, however, is unlikely 
because the surface brightnesses of the three components are significantly different from 
those observed in IRC+10216 and predictions of simulations. Moreover, the NE bow-shaped 
filament and arc have very faint counterparts on opposite sides of the central star (see the 
right panel of Fig. [TJ, suggesting that they were likely developed through different ejection 
events. The appeara nce of the NE a r c is sim ilar to that of the brightening shell as predicted 
by the simulations of IWareing et al.l (J2007aj), probably suggesting that it was formed by the 
interaction between a fast wind and the early AGB wind. Furthermore, the appearances 
of the NE arc and its SW counterpart suggest that it is li kely to be a debris of the ring, a 
feature that has been frequently detected in halos of PNs (ICorradi et al.l 120041 ) . 



E quations 1 and 2 are s teady-state solutions. In the simulations of Betelgeuse's bow 
shock, iMohamed et al.l ( 120121 ) found that it takes a long time (> 30 000 yr) for the bow shock 
to reach a steady state, before which the R(90°)/R(0°) ratio increases with time from 1 to 
y/3. Assuming that the symmetric axis of the halo of NGC 7293 lies in the sky plane (i.e. 
i = 0°), we obtain R(90°) / R(0°) = 1.32. Therefore, a possible explanation for the mismatch 
in the model and observation results is that the halo has not reached a steady state even 
though it is very old. 

In addition, Wilkin's solutions are based on the assumption of a constant stellar wind, 
which is applicable for the early AGB stage. However, as the star evolves up the AGB, 
both the mass- loss rate and wind velocity increase in magnitude. A jump in the stellar wind 
velocity in the post-AGB phase of evolution will cause interactions among stellar winds, 
which will also have an effect on the shape of the halo boundary. The shape of the halo can 
also be affected by additional mechan isms such as magnetic fields, which may distort the 
halo boundary. Ivan Marie et al.l (120 111 ) performed simulations including drag force between 
dust and gas. They found that depending on sizes of dust grains, infrared observations of 
bow shocks may not represent the morphologies of gas. The most likely explanation of the 
discrepancy between the model and the observed halo morphology is a latitude dependence of 
the density distribution of the AGB wind. AGB winds are known to develop non-spherically 
symmetric outflows near the end of AGB evolution, and this can easily change the shape of 
the interaction region. This equatorial enhanced mass loss have been attributed to binary 
central stars, star-planet interaction, star spots, and other mechanisms, but the exact cause 
is yet to be identified. Another possible complication is the angular dependence of the later- 
developed fast wind (bipolar outflow), which can distort the density distribution of the AGB 
envelope. Further numerical simulations including multiple stellar winds will be needed. 
Additional variables such as time- dependence and angle- dependence of outflows need to be 
investigated. The observed halo morphology presented here will serve as useful constraints 
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to these models. 



Recently, iRamos-Larios et al.l ( 120121 ) discovered an outer halo around the young PN 
IC418 based on the WISE images. Unlike that of NGC7293, the halo of IC418 does not 
exhibit a cometary shape and has a well defined boundary, suggesting that it has not been 
striped. Its elliptical a ppearance indicates that the density and velocity AGB winds are in- 
deed angle-dependent. IVillaver et al.l ( 120031 ) performed 2-D numerical simulations of the PN- 



ISM interaction considering variable AGB mass-loss rates. The shape of the NGC 7293 halo 



i s quali tatively similar to their results. The subsequent 3-D simulations by IWareing et al. 



( j2007at ) predicted four stages of the PN-ISM interaction. According to their classification 
scheme, this evolved PN is in t he first stage when the main nebula is unaffected by the inter- 
action. IWareing et al.l (j2007bl ) showed that the Kelvin-Helmholtz instabilities at the shock 



front can lead to vortices in the wakes of AGB stars. Although we do not see unambigu- 
ous vortex-like structures in the halo of NGC 7293, t he "SW arc" seemingly resembles the 
vortices reproduced b y the models of IWareing et al.l ( ]2007bl ). Through theoretical analysis, 
Dgani fc Sokerl (119981 ) found that because of the Rayleigh- Taylor instability the PN halo can 
be fragmen ted by the flow of ISM material . This phenomenon has been discovered in the PN 
NGC 6772 dRamos-Larios fc Phillip] bopgh . We note that the halo appearance of NGC 7293 
highly resembles to that of NGC 6772. Both show bright rimmed region in the shock head, 
fan-like material distribution in the tail, and ray-like enhancements in the interior. It follows 
that the Rayleigh- Taylor disruption is also present in the halo of NG C 7293. NGC 7293 is pre- 
suma bly more evolved than those PNs with halos investigated by IRamos-Larios &: Phillips 
( 120091 ). We do not find tight relations between the properties of the halos and the ages of 
their host PNs. This should be ascribed to different environments surrounding these PNs 
which play an important role in the PN-ISM interaction processes. 



5. CONCLUSIONS 

From the imaging data of the WISE all-sky survey we discover a mid-infrared halo with 
a diameter of ~ 40' around NGC 7293. The halo is totally invisible at optical wavelengths 
and emits light with a peak intensity near 12 /zm. The morphology of the halo suggests that it 
has been shaped by wind-ISM interaction, and the NE arc is brightened through wind-wind 
interaction. We find that the shape of the halo boundary cannot be explained by the analytic 



solutions of IWilkinl (119961 ) under the thin shock-shell and steady-state approximations. This 
suggests that mass loss during the progenitor AGB phase has a more complicated history 
due to time- and/or angular-dependent outflows. 

PN systems are the results of dynamical interactions between the halo (remnants of AGB 
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wind) with the interstellar medium, as well as the interactions between fast outflows with the 
remnants of the AGB wind (interacting winds). The quantitative information derived from 
this study is useful for the understanding of both types of dynamical interactions. NGC 7293, 
being the PN closest to us, is a valuable laboratory for the study of dynamic processes at 
work in the astrophysical environment. 

We thank Sze Ning Chong for helpful discussions. We also thank the anonymous referee 
for helpful comments. This paper makes use of data products from the Wide-field Infrared 
Survey Explorer, which is a joint project of the University of California, Los Angeles, and 
the Jet Propulsion Laboratory/ California Institute of Technology, funded by the National 
Aeronautics and Space Administration. Support for this work was provided by the Re- 
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Fig. 1. — Composite-color WISE image (left panel) and false-color GALEX image (right 
panel) of NGC7293. Left is east, and up is north. In the left panel, the WISE 3.4 /zm band 
is assigned the color blue, the 12 /im band the green, and the 22 /im ban d the red. The labels 
"NW outer feature" and "SE outer feature" are the terms assigned by lO'Dell et al.l (120041 ) 
for the bipolar outflow. In the right panel, the top-right corner panel shows a zoomed-in view 
of the black box. The cross and arrow at the center of each panel represent the location of 
the central star and the direction of proper motion. An axisymmetric halo is clearly visible 
in the WISE image. Note that in the GALEX image, the spike above the NW outer feature 
is artificial. 
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Fig. 2. — WISE 12 /im image shown in gray scale. The area of the image is 44' x 44' with 
the north direction 60° clockwise from the vertical axis. The cross and arrow represent the 
location of the central star and the direction of proper motion (as in Fig. [p. The boundary 
of the halo is indicated by the cyan dashed circle with a radius of 20' from its center (denoted 
by a cyan colored asterisk). The position of the NE arc is indicated by the green dashed 
circular arc with a radius of 11.7' from the central star. The red and blue dashed curves 
represent the four models listed in Table [Hand the number of the model is labeled next to 
the curve. The projected locations of the apexes of models 2 and 4 are marked by labels "a" 
and "b". 
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Table I. The Modeling Results". 



Model 


% 


Rq 


V*/y/niSM 




(°) 


(10 18 cm) 


(kms -1 cm 1 - 5 ) 


1 





2.78 


17.2 


2 


60 


1.80 


26.5 


3 





2.13 


22.3 


4 


60 


1.44 


33.0 



a Assuming M = lO _5 M yr^ 1 , v w = 
lOkms -1 , and the dimensionless mean 
molecular weight of 1.33. 



